Superconductor-Quantum Dot-Superconductor Junction in the Kondo Regime 
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Electron transport between two superconductors through an Anderson impurity in the Kondo 
regime is investigated within the slave boson mean field approximation. The current, shot noise 
power and Fano factor are displayed versus the applied bias voltage in the subgap region and found 
to be strongly dependent on the ratio between the Kondo temperature Tk and the superconducting 
gap A. In particular, the I — V curve exposes an excess current in the limit Tk/ A ^> 1. 



Background: A number of recently developed exper- 
imental techniques allow for detailed investigations of 
electronic transport through atomic-size metallic conduc- 
tors . Usually, transport properties of such systems are 
strongly affected by Coulomb interactions. Novel physi- 
cal effects emerge if electrodes of an atomic-size contact 
become superconducting. In that case the mechanism 
of multiple Andreev reflections |^ (MAR) plays a domi- 
nant role being responsible for both dc Josephson effect 
and for dissipative currents at subgap voltages. Further 
possibilities for experimental investigation of an interplay 
between MAR and Coulomb effects in systems with few 
conducting channels are provided by recently fabricated 
superconducting junctions with a weak link formed by a 
carbon nanotube [^j4|. 

Recently we developed a theory Q for the study of an 
SAS junction consisting of an interacting quantum dot 
{A) connected to superconducting [S) electrodes. It en- 
ables the analysis of MAR in superconducting contacts 
with few conducting channels in the presence of electron- 
electron interactions. An interplay between MAR and 
Coulomb effects is responsible for novel effects such as 
an interaction-induced shift of the subharmonic gap steps 
on the I — V curve and Coulomb blockade of MAR. The 
latter may result in a strong suppression of the subgap 
current through the dot at sufficiently low temperatures. 

The combined effect of MAR and electron-electron in- 
teractions on the shot noise |^ in superconducting quan- 
tum dots was studied in Ref. [Q. It was demonstrated 
that interaction effects can strongly suppress the shot 
noise power at subgap voltages. At the same time, the 
Fano factor (proportional to the ratio of the shot noise 
power and the current) was found to be nearly indepen- 
dent of interaction. 

Our previous analysis was restricted to physi- 

cal situations outside the Kondo regime implying that 
all relevant energies in the problem were taken much 
higher than the Kondo temperature Tk ~ 0. In the 
present work we study the opposite limit of sufficiently 
high Kondo temperatures. This case is relevant, e.g. 
for junctions composed of carbon nanotubes [^[Q] where 
the Kondo effect with rather high Kondo temperature 
Tk ~ 1.6 A' was recently observed M. Properties of 



superconducting quantum dots in this strong coupling 
regime are entirely distinct from those exposed earlier 
in the limit Tk 0. On a qualitative level, this dif- 
ference between weak and strong coupling regimes can 
be briefly summarized as follows: In the weak coupling 
limit (low Tk), although the number of Andreev reflec- 
tions n 2A/eV may be large at low voltages V, both 
the current and shot noise power remain rather weak. 
This is due to the low effective transparency F of the 
junction as a consequence of strong repulsive electron- 
electron interaction (Coulomb blockade). Large-n pro- 
cesses are therefore damped as F". By contrast, in the 
strong coupling Kondo regime the effective transmission 
is much larger and a ballistic-like channel opens up in- 
side the dot. Hence, an interplay between MAR and the 
Kondo resonance is expected to yield an excess current 
in the I ~V curve, similarly to the case of noninteracting 
ballistic junctions. At very large values of Tk and in the 
low voltage limit this current should approach the nonin- 
teracting result Q Jar = 4eA//i. Analogously, the shot 
noise power is expected to display a pronounced maxi- 
mum a.t V — and should decay as 1/V at small bias as 
is familiar in the standard noninteracting SNS junction 
|]l0| . Below we will present a quantitative analysis which 
fully supports this qualitative physical picture. 

Model and effective action: The pertinent system is 
represented by two half planar electrodes {L and R) sep- 
arated by the line a; = 0, and weakly coupled to a point- 
like Anderson impurity A located at the origin. This 
model is of interest, for instance, in connection with re- 
cent experiments [pd],r2| on semiconductor quantum dots. 
It was shown there that tunneling takes place through a 
separate state with features of a Kondo behavior (a tun- 
able Kondo effect). 

The system dynamics is governed by the Hamiltonian 



H = Ht 



Hd + Hf + Hr 



(1) 



in which i?L._R are the BCS Hamiltonians of the elec- 
trodes which depend on the electron field operators 
'0L(/s)tT(r, f) where r = {x,y) and cr = ± is the spin 
index. As in Refs. the dot is described as a sin- 
gle level Anderson impurity A with energy eg < and 
Hubbard repulsion parameter U. In the Kondo regime 
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of interest here we set C/ — > cxd and assume |eo| to exceed 
any other energy scale except U. In this case it is conve- 
nient to express the dot and the tunneling Hamiltonians 
Hd and Ht via slave boson (operators 6,6^) and slave 
fcrmion (operators c,c^) auxiliary fields JTst . Explicitly, 
Hd = eoE^4c<T and Ht = T X^^^ ct^V'jTlO, i)+h.c., 
where T is the tunneling amplitude. Finally, the Hamil- 
tonian of the system must also include a term which pre- 
vents double occupancy in the limit [/ — > oo. This term 
reads He — '4r'^cr + b^b — 1), where A is a Lagrange 

multiplier. 

Let us now consider the dynamical "partition function" 



V[F] exp(i5). 



(2) 



where the path integral is carried out over all fields [F] 
and the action S is obtained by integrating the La- 
grangian pertaining to the Hamiltonian (^ along the 
Keldysh contour. The procedure of integrating out the 
electron fields of the bulk electrodes V'L(fl)o-(r, was de- 
scribed in details in ||^ . As a result we arrive at the effec- 
tive action expressed in terms of the Green functions of 
the bulk superconductors. Our next step is to integrate 
out the variables corresponding to the Fermi operators cj. 
and Co- of the dot. The corresponding integral is Gaus- 
sian, which yields, 



5eff = -iTi-lnG-^ - / dt[Xa^{bb - I)]. 



(3) 



Here A = (Ai,A2), b = (61,62) and az are diagonal ma- 
trices acting in Keldysh space. Similarly to Ref. ||^, the 
inverse propagator G^^ depends on the Green functions 
of the electrodes. 

Mean field slave boson approximation (MFSBA): In or- 
der to describe the Kondo regime we will treat the slave 
boson fields 61 and 62 in Eq.(^) within the dynamical 
mean field approximation. Performing the variation of 
the effective action with respect to 6i_2 and Ai^2 and then 
setting 61 = 62 = 6 and Ai=A2 = A we arrive at two self- 
consistency equations that determine the parameters 6 
and A. 

Before presenting these equations let us specify the ex- 
pression for the inverse propagator . Performing the 
standard basis rotation in Keldysh space one finds 

r62 

G-i(6, e') = 5{e - e'){e ~ r,e) + — r,g+(e, 6')t„ (4) 

where e = eo + is the renormalized level position (in 
the Kondo limit one has e ~ 0) and F oc is the 
usual transparency parameter. Here and below we de- 
fine g± = QL^gR, where 



' 27r 



(5) 



are Keldysh matrix Green functions of left and right elec- 
trodes and 'p/e = V{t) is the bias voltage across the dot. 



The matrix g has the standard structure with retarded 
and advanced Green functions 



(e±zO) + |A|T^ 
V(e±*0)2-|A|2' 



(6) 



as diagonal elements g^/^ and the Keldysh function 
g^(e) = (g«(e) - g^{e))ta.Dh{e/2T) as the only nonzero 
(upper) off-diagonal element. The Pauli matrices Tx.y,z 
act in Nambu space. The inverse of the matrix (^) is 
formally performed, leading to a 2 x 2 Keldysh Green 
function with three elements, 



G 



R.A 



-Tzgr Tz\ , 
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(7) 
(8) 



In order to explicitly write down the self-consistency 
equations let us introduce the bare Kondo temperature 
= i:)ea;p[-7r||eo|/(2F)] and define T9 = T^X, where 
D is the energy bandwidth. Then our MFSBA equations 
take the form 



X = -^TrG^r,, 



2t: 



K 



A = ^Tr[G^'(5« + 5^) + (G^ + 0^)9""' 



(9) 
(10) 



where the trace also includes energy integration. Eq. (||) 
effectively determines the Kondo temperature (through 
the parameter X\ and reflects the constraint which pre- 
vents double occupancy in the limit U 00. The second 
self-consistency equation (^0|) defines the renormalized 
energy level position e. 

Let us briefly discuss the validity range of the present 
analysis. The MFSBA is known to encode the Kondo 
Fermi-liquid behavior at low temperatures. An impor- 
tant parameter here is the ratio between the Kondo 
temperature and the superconducting gap tj^ = T^/A 
|p^-16|. For i/f > 1 a Fermi liquid behavior is expected. 
Accordingly, in this regime Eq. (^ should have a nonzero 
solution X which corresponds to nonzero Tk- On the 
other hand, in the limit of large A the only possible so- 
lution is a trivial one 6 = (and, hence, Tk = 0) In 
this case - as it was demonstrated in Ref. Q - the prob- 
lem can be treated within the dynamical mean field ap- 
proximation for the bare Anderson Hamiltonian. Quanti- 
tatively, the MFSBA is reliable only for sufficiently large 
values of tx- We believe, however, that it can provide 
useful qualitative information also for moderate values 
of tfc describing a crossover between the Kondo regime 
and the Coulomb blockade behavior It is worth not- 
ing here that the applied bias voltage V also attenuates 
the Kondo resonance and lowers Tk- Hence, for the relia- 
bility of the MFSBA in non-equilibrium situations, both 
A and eV should not exceed the Kondo temperature. 
Attention below is mainly focused on the subgap voltage 
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regime eV < A in which case Ik appears to be the only 
relevant parameter. 

/- V curve: The standard expression jlTt for the tunnel- 
ing current operator between the dot and one (e.g. the 
right) electrode reads, 



(1:2) 



= ±- 



■En 



l±cr. 



V'i?.fe(0)-h.c.]. 



(11) 



Here ipRkiO) is the Fourier transform of '0jj(O, y, t) with 
respect to y and (1,2) refer to Keldysh indices. As be- 
fore, it is convenient to integrate out the ip-Helds and 
express the current / through the dot in terms of the 
Green functions of the bulk electrodes. This procedure 
has been described in details in Ref. As a result we 
obtain 



(12) 



where we defined g — glr^^ — r^gC: and g ' ~ g 



~R,A 



-R,A 



Tzg'S"- Being combined with eqs. (^, (10) the result 
( p^ ) can be conveniently used for computing the trans- 
port current of an SAS junction in the Kondo regime for 
different values of Ik- 

For sufficiently large Ik we anticipate a strong Kondo 
resonance and the I — V curve is expected to resemble 
that of purely ballistic junctions without interaction [Q. 
Indeed, in the limit of large 3> 1 which corresponds 
to the unitary, pure ballistic case, our expression for the 
current (|l2|) reduces to that of Ref. [D . This agreement is 
further supported by our numerical analysis carried out 
for tK — 100 (dashed-dotted curve in Fig.l). Calculation 
of the current was also performed for tx = 5 and 1.6 and 
the results are displayed in Fig. 1. 




eV/A 



curves correspond respectively to tx = 100, 5 and 1.6. 



For tK = 5, a, pronounced excess MAR current is 
clearly exposed in the I — V curve, though its magni- 
tude turns out to be smaller than in the unitary limit 
Ik ^ 1- One also observes sub- harmonic MAR steps 
(which are hardly visible in the case Ik ^ 1). For a 
lower value tj^ = 1-6 the current is strongly suppressed 
(dotted line in Fig. 1), and the I — V curve resembles 
that of a low transparency SAS. MAR steps in the I — V 
curve become more pronounced as compared to the case 
of higher tK ■ 

Let us briefly summarize our results for the current. In 
the limit of large tK the I — V curve is practically inde- 
pendent of tK and resembles that of a ballistic junction, 
as indicated in the upper curve in Fig.l. This effect is 
physically similar to that discussed in Refs. |l9|j2C|l where 
stimulation of the dc Josephson current was found in the 
unitary limit of the Kondo regime. For lower values of 
tK the junction is driven away from the unitary limit and 
the system behavior becomes much richer. It reflects the 
influence of both A and V on the Kondo resonance and 
on the actual value of the Kondo temperature. While for 
tK = 5 some imprint of ballistic behavior still remains 
(finite excess current), the I — V curve noticeably devi- 
ates from that obtained in a non-interacting limit. For 
tK = 1.6 the competition between gap-related suppres- 
sion of the Kondo effectand the effective transparency of 
the junction becomes essential leading to further decrease 
of the current. For even lower values tK, the Kondo 
physics is no more relevant and one crosses over to the 
Coulomb blockade regime 

Shot noise: The shot noise spectrum is usually defined 
as the symmetrized current-current correlation function 
1^ . Being expressed via the operators ( pT| ) it reads 



K{t,,t2) ^ h[{fi('Ht,)i(^Ht2)) 

+ {fl^'\t2)I^^Hti))-2{I)% 



(13) 



where T is the time ordering operator and (...) denotes 
quantum averaging with the Hamiltonian (|^) . Substitut- 
ing = (4^'^^ - ll^;^^)/2 into eq. (^ we obtain an 
expression for K{ti,t2) which involves integration over 
surface fields and dot electron slave particle field opera- 
tors. The first integration involving the Green function 
matrix g is Gaussian and can be done exactly. Integra- 
tions over the dot slave fermion fields is completed within 
the dynamic MFSBA. After Fourier transform with re- 
spect to ti — t2 it is possible to express the power noise 
spectrum K(lu) in terms of the Green functions of the 
entire system (0), (^. It is convenient to decompose 
K = {Ki + K2}?A/{8h) with the result. 



Xt 



TT{{gf gf){G'' G^) - G^}, (14) 



Fig. 1 The averaged current / (in units of AeA/h) versus 
the bias V (in units of A/e) for an SAS junction at sub-gap 
voltages with F/Tj^ = 200. Dashed-dotted, solid and dotted j^, 



(XtK) 



Tr{(G^.g)2 - 2r,gr,G^gG« - 
o 

[2gG«r,g5G^' + (ff^G^)^ - {G^g'^Tzf + h.c.]}. (15) 
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Expressions (supplemented by the self- 

consistency eqs. (I), (p^) are then solved numerically 
for the same set of parameters F/T]^ = 200, = 100, 5 
and 1.6. The results for the shot noise power spectrum 
K versus the applied voltage V are displayed in Fig. 2. 
These results are clearly correlated with those for the 
I — V curve and can be summarized as follows. 

In the limit i/f ^ 1 our results are consistent with 
those obtained for purely ballistic junctions ||l^. In par- 
ticular, we mention that the noise spectrum K exhibits a 
1/V dependence at low voltages. At lower Ik the physics 
is distinct. For tx — ^ (solid curve) the noise spectrum 
still shows features typical for a junction with relatively 
high transparency (see |^ and Fig. 1 therein), while 
the results for Ik — 1.6 (dotted curve) are more similar 
to those for a low transparency junction. For tx = 5 
we get a pronounced noise power at low voltage whereas 
at higher voltages eV ^ A the shot noise is significantly 
suppressed. At low voltage K scales approximately as 
l/V. The sequential tunneling picture is not valid in the 
unitary limit ti<- 3> 1 as well as for the intermediate value 
tK — 5 due to the interference between different MAR 
processes. Yet, this picture is partially restored at lower 
Kondo temperatures, as is demonstrated by our results 
obtained for tx = 1.6. 



1.0 
eV/A 



Fig. 2 The shot noise power K (in units of 2e^A/h) as a 
function of V (in units of A/e) for an SAS junction. The 
parameters and notations are the same as in Fig. 1. 
This effect of sequential tunneling due to MAR is most 
clearly manifested through the Fano factor K/2eI de- 
picted in Fig. 3 as a function V. 
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Fig. 3 The Fano factor K/2eI versus V. The parameters 
and notations are the same as in Fig. 1. 

We observe that for tx = 5 the Fano factor at eV > A 
is substantially lower than the expected value K/2eI — 2 
originated from Andreev reflections. On the other hand, 
for tx — l-Q the Fano factor is closer to this value. 

In conclusion, we have analyzed an important physical 
problem involving strong correlations, the Kondo effect 
and superconductivity. These aspects can be combined 
in an SAS junction consisting of an Anderson impurity 
(in the Kondo regime) located between two supercon- 
ducting electrodes, which is experimentally feasible. We 
have developed a theoretical framework by which it is 
possible to investigate an interplay between MAR and 
Coulomb effects in the Kondo regime T < A < Tk- In 
this limit we have exposed the nonlinear I—V characteris- 
tics and calculated the shot noise power spectrum of SAS 
junctions at subgap voltages eV < 2A. We have found 
that at sufficiently large the Kondo resonance plays 
the dominant role effectively making the junction behav- 
ior similar to that of highly transparent non-interacting 
weak links . This physical situation is qualitatively 
different from the Coulomb blockade regime encountered 
in the limit A > T > Tk which we have analyzed in 
our previous works [^,0 . A crossover between these two 
physically different regimes occurs at Tk ~ A and is also 
- at least qualitatively - described within our theoretical 
framework. 
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